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ABSTRACT 

We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to 
map the gas kinematics of the inner ~ 200x500pc of the Seyfert galaxy NGC 4151 
in the Z, J, H and K bands at a resolving power ^ 5000 and spatial resolution of 
~8pc. The ionised gas emission is most extended along the known ionisation bi- 
cone at position angle PA=60-240°, but is observed also along its equatorial plane. 
This indicates that the AGN ionizes gas beyond the borders of the bi-cone, within 
a sphere with w 1" radius around the nucleus. The ionised gas has three kinematic 
components: (1) one observed at the systemic velocity and interpreted as originating 
in the galaxy disk; (2) one outflowing along the bi-cone, with line-of-sight velocities 
between —600 and 600 kms -1 and strongest emission at ± (100 — 300) kms" 1 ; (3) and 
another component due to the interaction of the radio jet with ambient gas. The radio 
jet (at PA=75-255°) is not aligned with the NLR, and produces flux enhancements 
mostly observed at the systemic velocity, suggesting that the jet is launched close to 
the plane of the galaxy (~ plane of the sky). The mass outflow rate, estimated to 
be «lM Q yr _1 along each cone, exceeds the inferred black hole accretion rate by a 
factor of ~ 100. This can be understood if the Narrow-Line-Region (NLR) is formed 
mostly by entrained gas from the circumnuclear interstellar medium by an outflow 
probably originating in the accretion disk. This flow represents feedback from the 
AGN, estimated to release a kinetic power of E w 2.4 x 10 41 ergs _1 , which is only 
~ 0.3% of the bolometric luminosity of the AGN. 

There is no evidence in our data for the gradual acceleration followed by gradual 
deceleration proposed by previous modelling of the [OIII] emitting gas. Our data allow 
the possibility that the NLR clouds are accelerated close to the nucleus (within O'.'l 
- «6pc) after which the flow moves at essentially constant velocity (« 600 km s -1 ), 
being consistent with NIR emission arising predominantly from the interaction of the 
outflow with gas in the galactic disk. 

The molecular gas exhiibits distinct kinematics relative to the ionised gas. Its 
emission arises in extended regions approximately perpendicular to the axis of the 
bi-cone and along the axis of the galaxy's stellar bar, avoiding the innermost ionised 
regions. It does not show an outflowing component, being observed only at velocities 
very close to systemic, and is thus consistent with an origin in the galaxy plane. This 
hot molecular gas may only be the tracer of a larger reservoir of colder gas which 
represents the AGN feeding. 

Key words: Galaxies: nuclei, Galaxies: active, Galaxies: Seyfert, Galaxies: kinemat- 
ics, Galaxies: jets, Galaxies: individual (NGC 4151) 
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1 INTRODUCTION 

This is a continuing study of the narrow-line region (here- 
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tained with the Gemini Near-infrared Integral Field Spec- 
trograph (NIFS). The data comprise spectra of the inner w 
200 x 500 pc 2 , at a spatial resolution of ~ 8 pc at the 
galaxy, covering the wavelength range 0.95-2.51 fim at a 
spectral resolving power R ^ 5000. In a previous paper 
l|Storchi-Bergmann et al.l 120091 . hereafter Paper I) we have 
used these data to map the NLR intensity distributions 
of 14 emission lines, as well as their ratios. The main re- 
sults were the distinct flux distributions and physical prop- 
erties observed for the ionised, molecular and coronal gas. 
The ionised gas is co-spatial with the ionisation bi-cone ob- 
served in the optical [Qui] emission at position angle (here- 
after PA) 60° (|Evans et al.lll993l ; lHutchings et al.lll998l), and 
seem s to trace the outflow along the bi-cqne (|Das et al.l 
l2005l ; ICrenshaw ei~al1 12000| ; iHutchings et all Il999h . In the 
inner region, the NIR ionised-gas emission extends beyond 
the borders of th e cone, which does no t have a sharp apex 
as noted also by iKraemer et all (120081 ). The H 2 molecular 
gas intensity distribution, on the other hand, avoids the re- 
gion of the bi-cone and seems to originate in the galaxy 
disc, while the coron al gas emission is b arely resolved. In 
another recent work (|Riffel et al. 2009b;), we have studied 
the unresolved nuclear continuum showing that its origin 
is emission by hot dust within ~4pc from the nucleus - 
as expe cted from the dusty torus postulated by the unified 
model l|Antonucci fc Millerlll985h . The origin of this struc- 
ture is probably a dusty wind that originates in the outer 
parts of the accretion disk. Thi s wind is probably clumpy , 
as suggested by recent models (lElitzur fc Shlosmanl 12006). 
and necessary in order to allow the escape of radiation along 
the equatorial plane of the bi-cone in order to ionize the gas 
and produce the observed flux maps of Paper I. 



In the present paper we use the NIFS data to map the 
NLR kinematics of NGC4151. Although many papers have 
been devoted to such a study, most of them are based on 
long-slit spectroscopy ob tained in the optical with the Hub- 
ble S p ace Telescope (e . g.lWinge et al.|[T9 97: Hutc hings et al.l 
1 19981 ; ICrenshaw et al.ll2000l ; iDas et al.ll2005l ). The present 
study was performed in the near-IR, a waveband which is 
less affected by dust, and using integral field spectroscopy, 
which allows a full two-dimensional (hereafter 2D) coverage 
of the NLR kinematics. With our data and analysis we aim 
to examine previous claims of acceleration of the gas along 
the NLR, quantify the mass outflow rate and the correspond- 
ing feedback power, as well as investigate the origin of the 
NLR gas. As NGC4151 harbors the closest bright AGN, its 
NLR is one of the best suited for this type of study. 



Our approach for the analysis of the NLR kinematics of 
NGC4151 is as follows. First we map the centroid velocity 
and velocity dispersion of the ionised, molecular and coronal 
gas, obtained from fits to the emission-line profiles and com- 
pare our results with those of previous studies. Then we use 
an alternative approach to map the NLR kinematics, made 
possible by integral field spectrographs: a "velocity tomog- 
raphy" of the emitting gas, obtained by slicing the line pro- 
files in velocity bins of 60kms _1 , producing channel maps 
which provide a clearer view of the velocity distribution in 
the different gas phases. We have grouped the channel maps 



together in a sequence of velocity bins, generating movies 
which can be recovered at the authors' website. 

We adopt in this paper the same distance to NGC 4151 
as in Paper I: 13.3 Mpc, corresponding to a scale at the 
galaxy of 65 pc arcsec -1 (|Mundell et al.ll2003h . The paper is 
organised as follows. In §[5] we provide a quick description of 
the data, in §[3]we present the centroid velocity and velocity 
dispersion measurements, in §[4] we present the emission- line 
"tomography" , in §[S] we discuss our results, in §[B] we es- 
timate the mass outflow rate along the NLR and compare 
with the mass accretion rate and in §[7] we present our con- 
clusions. 



2 THE DATA 

The observations comprise three pointings obtained with 
the Gemini North te lescope with the NIFS instrument 
|McGregor et all [20031 ) operating with the ALTAIR adap- 
tive optics module, each covering a square field of view 
of ~ 3.0 x 3.0arcsec 2 . One of the NIFS fields was cen- 
tred on the NGC 4151 nucleus and the others were offset 
by ± 2.5 arcsec along position angle 75°, resulting in a total 
field of view of 8.0 x 3.0 arcsec 2 , with the longest extent ori- 
ented approximately along the radio jet (I Pedlar et al.|[l993l ; 
iMundell et"aT]|l995l . liool ). 

The spectra cover the Z, J, H, and K spectral bands 
at two-pixel resolving powers of 4990, 6040, 5290, and 5290, 
respectively. The angular resolution, as determined from the 
FWHM of the PSF, ranges from 0.12 ± 0.02 arcsec in the H 
and K spectral bands to 0.16 ± 0.02 arcsec in the Z and J 
bands. This corresponds to a spatial resolution at the galaxy 
of 8 ± 1.3 pc at H and K bands, and 10 ± 1.3 pc in the Z and 
J bands. 

More details about the data and the data reduction 
can be found in Paper I, as well as flux measurements of 61 
emission lines from three characteristic locations along the 
NLR. In the present paper, we concentrate on the kinematics 
measured for the strongest emission lines of the ionised and 
molecular gas, namely: [S in] A 0.9533 /im, He I A 1.0833 /im, 
[P n] A 1.1886 /im, Pa/3 A 1.2822 ^m, [Fe n] A 1.6440 /im and 
H2 A 2.1218 /J,m. We have also measured the kinemat- 
ics for the coronal lines of [S ix] A 1.2523 /im and 
[Cavin] A 2.3220 /im. We first present and discuss the cen- 
troid velocity and velocity dispersion maps and then present 
the velocity tomography. 



3 CENTROID VELOCITY AND VELOCITY 
DISPERSION MAPS 

Centroid velocity and velocity dispersion (a) maps were 
obtained by fitting Gaussians to the emission-line profiles. 
The velocities were obtained from the central wavelength 
of the Gaussians (after subtraction of the systemic veloc- 
ity of 997kms _1 , and applying a heliocentric correction of 
+ 17kms _1 ) while the velocity dispersions were obtained 
from the full width at half maximum (hereafter FWHM) 
of the Gaussians as a =FWHM/2.355. 

In many locations along the NLR of NGC 4151, the 

1 http:/ /www. if.ufrgs.br/~thaisa/ifu_movics/ngc4151 
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strongest ionised gas emission lines - those of [S ill] and He I 
- show two kinematic components, with a mean velocity sep- 
aration between them ranging from 200 to 350 km s" 1 . We 
were successful in fitting two Gaussian components to both 
these lines, although the He I line has a broad component 
and a nuclear absorption component which introduces un- 
certainties in the fitting procedure for regions close to the 
nucleus. In a few regions, we observe double components also 
in [Fell], Pa/3, and Br7. Elsewhere, these lines are broader 
at locations where double components are observed in [S ill] 
and He I. This suggests that the [Fell], Pa/3, and Br7 emis- 
sion lines are also double. The lower signal-to-noise ratio 
(hereafter S/N) of these lines, the presence of a broad com- 
ponent for the H I emission-lines close to the nucleus, and the 
fact that the [Fen] emission line appears to be intrinsically 
broader contribute to blend together two distinct kinematic 
components. 

Double components have been previou sly observed i n 

the 
and 



O nil A 5007 A emissio n -line profil e s bvlSchulj l|l990l ). 
Crenshaw et all (|2000T ). iDas et al.1 (120051 ) using higher 



resolution (R~ 9000) spectra, showed that in some loca- 
tions there are even three or four components. Our lower 
spectral resolution did not allow us to resolve more than 
two components. 

Typical fits to the [S ill] emission-line profile are shown 
in Fig. [1] corresponding to locations (A, B, C, D and E) 
identified in the [S ill] radial velocity maps which are shown 
in Fig. Because the S/N is not high enough to allow the 
fit of each component with the necessary number of free pa- 
rameters, we had to constrain the width of the two Gaussian 
components to the same value. The implication of this as- 
sumption is that there will be some uncertainty in the gas 
kinematics derived in this way, mainly in the velocity dis- 
persions, as discussed bellow. Although in the next section 
we derive the kinematics in terms of channel maps, with no 
assumptions, we decided to include the kinematics in terms 
of two components in order to be able to compare our re- 
sults with those of previous studies, in which this was the 
method used to derive the gas kinematics. In addition, the 
separation in two components helps in the identification of 
the origin of the gas emission. 

In the top panels of Fig. [2] we show velocity maps in 
the emission lines of [S ill] and He I of regions with one kine- 
matic component, after subtraction of the systemic velocity. 
In the middle and bottom panels we show the velocity maps 
of the regions with two kinematic components: the middle 
panel shows the high velocity component, corresponding to 
negative (blueshifted) velocities to the SW and positive (red- 
shifted) velocities to the NE, while the bottom panel shows 
the low velocity component. 

For the [S ill] emission, the lowest velocities are observed 
along the galaxy minor axis (PA ~ —68°) in the upper panel 
of Fig. [2] (single component), and in neighboring locations 
where two components were fit to the data. This low velocity 
component, essentially at zero velocity, probably originates 
in the galaxy disc, while the redshifted and blueshifted ve- 
locity components originate in the outflow along the bi-cone. 
The highest blueshifts are observed in the bottom part of the 
strip labelled "C" in the middle left panel of Fig. [2] while the 
highest redshifts are observed in the region labelled "B" . 

In order to better illustrate the magnitude of the ve- 
locity variations, and to investigate further the origin of the 



different gas components, we show in Fig. [3] the velocities 
as a function of distance from the nucleus along PA= 60°. 
The black dots represent the single components and the red 
dots represent the double components, while the labels al- 
low the comparison with Fig. [2] Fig. [3] suggests the presence 
of a ~ zero velocity component observed between ~ — 1'.'5 
and « l'.'5 from the nucleus along PA« 60°, whose origin 
is probably the galaxy plane, as pointed out above. Then 
there is a sequence of blueshifts at « — 250kms _1 to the 
SW (positive distances in Fig[3]) which probably originates 
in the far side of the approaching cone (see Fig. [9] and dis- 
cussion in Sec l5.3p . A similar sequence of redshifts at some- 
what lower velocities, up to ~ 150kms _1 to the NE can be 
identified with the near side of the receding cone. The high- 
est blueshifts to the SW and highest redshifts to the NE 
can then be identified with the near side of the approach- 
ing cone and the far side of the receding cone, respectively. 
These highest velocities thus originate in the walls of the 
bi-cone that are most inclined relative to the line-of-sight. 

Velocity dispersion (a) maps of the [S ill] and He I emis- 
sion lines for the single and two-component fits are shown in 
Fig. [4] We note again that the measured velocity dispersions 
are upper limits to the intrinisic gas velocity dispersion in 
regions where a single Gaussian component has been fit, pos- 
sibly to multiple unresolved components. For example, in re- 
gion E in Fig. 2] the velocity dispersion is ^ 150kms~ . But 
this large velocity dispersion is probably due to the blend- 
ing of at least two components, as it surrounds a region with 
two components (region D). Similarly, the measured values 
are only indicative of the gas velocity dispersion in regions 
where two components have been fit, as their widths have 
been constrained to the same value. This is the case of re- 
gion E, where there are a component at zero velocity and 
a blueshifted component, which do not necessary have the 
same velocity dispersion. Nevertheless, these a values should 
not be too far off because they are of the order of the a value 
of the zero velocity single component along the minor axis 
of the galaxy. 

The lowest a values for the single components are ~ 
100 km s" 1 . These are seen mostly along the minor axis of 
the galaxy. The low velocity dispersion suggests that the 
emitting gas in these regions is not disturbed and proba- 
bly lies close to the galaxy plane. The velocity dispersion 
is larger in other regions. This could be due to turbulence 
caused by the outflow or to multiple velocity components. 
The a values for the double components are in some regions 
(e.g. region "D") almost as low as the lowest values of the 
single components, suggesting that at least one of the com- 
ponents originates in the galaxy plane. Regions of higher 
a values in the two components are those labelled "B" - 
in the NE outflow - and the bottom part of the strip "C" , 
where the highest blueshifts are seen in the centroid velocity 
maps, supporting an origin in the front side of the SW cone, 
as discussed above. 

The centroid velocity distribution of H2 emission lines 
is dominated by velocities close to systemic. This emission is 
concentrated along the minor axis of the galaxy, whose orien- 
tation approximately coincides with that of the galaxy bar. 
There are some redshifts to the NE mainly along the galaxy 
major axis and some blueshifts to the SW along both the 
galaxy major axis and the cone axis (see Fig. [8} . The velocity 
dispersions are mostly constant at «80 kms -1 supporting 
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Figure 1. Typical profiles and fits to the [S III] emission line for the regions identified by letters A— E in Fig. [2] 



an origin for this gas in the galaxy plane. The coronal line 
flux distributions are just barely resolved and their centroid 
velocities are also close to systemic, while their velocity dis- 
persions are « 100 kms -1 . 



4 EMISSION-LINE SCANNING 

We have also mapped the gas kinematics by integrating the 
fluxes within wavelength slices extracted along the emission- 
line profiles corresponding to velocity bins of 63 kms -1 . 
These "channel maps" display the distributions of the gas 
emission at different velocities, which occur at different 
depths along the line of sight. In this sense, they can 
be described as velocities "tomography" of the NLR us- 
ing emission-line scanning. The channel maps are shown in 
Figs. [5] |6] [7] and [8] for the emission lines of [S m] , [Fell], Br7 
and Eb, respectively. Channel maps in the other ionised gas 
emission lines, such as He I and Pa/3, are similar to those in 
[S ill] and thus are not shown here. The lowest flux shown in 
the channel maps is 3 times the root mean square variation 
of the surroun ding sky flux. The bi-cone axis at PA= 60° as 
determined bv lDas et al.l (|2005l ) is represented by a dashed 
line in the figures, while the major axis of the galaxy at 
PA= 22° is represented by a continuous line. The contours 
of a radio MERLIN image (kindly made available by C. 
Mundell) have been overplotted on the central velocity bins 
(velocities close to systemic). 

In the channel maps of the [S ill] emission line, shown 
in Fig. [5] the most blueshifted gas at ~ — 600 kms -1 lies 
very close to the nucleus. For blueshifts of » —400 kms -1 
the emission extends further out, and for lower velocities 
broadens in the shape of a projected cone extending from 
the bi-cone axis towards the galaxy major axis. At the slice 
corresponding to —221 kms -1 there seems to be a spatial as- 
sociation between the line emission and the radio structure. 



In particular, both of them broaden towards the major axis. 
At velocities close to zero, and up to redshifts « 50 kms -1 , 
the emission is extended in all directions around the nucleus, 
evidencing the absence of a sharp cone apex. The largest 
elongation of this emission occurs along the horizontal axis, 
which is approximately the orientation of the radio axis. The 
overplotted radio contours show that, although there is not 
a detailed correspondence between the [S in] and radio im- 
age, the brightest regions of the [S ill] image coincide with 
the two inner knots of the radio image. For redshifts larger 
than 100 kms -1 , the gas is observed mostly to the NE in 
the opposite direction to that of the blueshifted gas, show- 
ing also a similar flux distribution: for the lowest redshifts 
the emission is distributed in a fan structure extending from 
the major axis of the galaxy to the bi-cone axis, while for 
the intermediate and the highest redshifts the emission is 
more elongated and concentrated closer to the bi-cone axis 
to the NE. 

Channel maps for the [Fell] emission line are shown 
in Fig. HJ] In the high-velocity blueshifted channels, this 
emission seems to be more collimated along the bi-cone 
axis than for other emission lines. Some blueshifts from 
« -221 kms -1 to zero velocity are observed also to the 
NE, unlike the slices in the other emission lines. The flux 
distributions for velocities close to zero and redshifted gas 
are similar to those in [Sin], although the highest redshifts 
are as high as the highest blueshifts and are observed up to 
similar distances from the nucleus. In all velocity channels, 
there seems to be more structure in [Fell] than in [Sill], in 
the form of knots of emission in closer association with the 
radio knots. This is particularly clear at the velocity bins 
from ~ —221 to 94 kms -1 : at —221 kms -1 the strongest 
emission traces the "flaring" of the radio structure to SW 
towards the galaxy major axis; between velocities of —45 and 
+45 kms -1 , the strongest emission approximately traces the 
two central radio knots and, for higher redshifts, there is an 
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Figure 2. Radial velocity maps for the single (top), and double components (middle and bottom panels) for [S III] (left) and He I (right). 
In the middle panel we show the high-velocity component, and in the bottom the low-velocity component. The continuous line shows 
the orientation of the galaxy major axis and the dashed line shows the orientation of the bi-cone. Emission-line profiles and fits typical 
of the regions marked with letters A— E are shown in Fig. [T] for the [S III] emission line. 



enhancement of the [Fell] emission at ~ 1 arcsec NE, at the 
location of another radio knot. 

Fig.[7]shows the channel maps in the Br7 emission line. 
These show similar morphologies to those seen in [S ill] , but 
with more structure. This can be attributed, at least in part, 
to the better spatial resolution achieved by the adaptive 
optics system in the K band relative to the Z band. 

The channel maps in the H2 emission line, shown in 
Fig. [8] are completely different from those of the ionised 
gas, in support of a di stinct or igin for the molecular gas, as 
already pointed out in |Paper~H . Most of the emission occurs 
in the ~ zero velocity channels. A small excess of blueshifts 
to the SW (from ~ 100 km s -1 ) and similar redshifts to the 
NE suggest rotation around the minor axis. The flux distri- 
butions show a hollow core around the nucleus, and beyond 
this region the emission surrounds the minor axis of the 



galaxy in the shape of two arcs extending up to 1 arcsec to 
NW and SE of the nucleus. A comparison with the channel 
maps in the ionised gas at the same velocity bins (around 
zero velocity) shows that the molecular gas emission avoids 
the innermost region of ionised gas, and the hollow struc- 
ture can thus be explained by the presence of ionised gas 
surrounding the nucleus at zero velocity. 



5 DISCUSSION 

5.1 Velocity variations as a function of distance 
from the nucleus 

In the centroid velocity maps for the ionised gas emission 
(Fig. El, the lowest velocities are observed at the nucleus, 
while the highest are observed away from the nucleus. This 
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is the case for NGC 4151 and also for 6 other Seyfert galaxies 
for which we have recently obtained kinematic data with 
the IFU of the Gemin i Multi-objetct Spectrograph (CMOS) 
|Barbosa et al.ll2009l ). 

The channel maps, on the other hand, show a differ- 
ent behavior : the highest velocit ies are usually observed at 
the nucleus l|Barbosa et alj|2009h . In the case of NGC 4151, 
blueshifts of w -600 km s" 1 in [S in] and [Fe n] are observed 
at the nucleus and in a narrow region along the bi-cone axis 
extending from the nucleus almost to the outer edge of the 
NLR (Figs. [5] and |6]) . Lower velocity blueshifted gas appears 
to fill the bi-cone and extend into the region of the radio 
jet near the systemic velocity. A similar structure is seen in 
redshifted gas with the emission at moderate redshifts ap- 
pearing to fill the receding bi-cone. Emission at the highest 
redshifts of > 400 kms -1 occurs at the nucleus and in a lo- 
calised region ~ Of! 9 to the NE of the nucleus that is offset 
from the bi-cone axis. This morphology differs from that of 
the highest velocity blueshifted gas and may simply reflect 
the clumpy nature of the NLR emission. 

We attribute the difference between the kinematics 
probed by the centroid velocity and channel maps to the 
fact that the centroid velocity probes the brightest emission, 



while the channel maps also probe fainter emission (in the 
wings of the line profiles). In the vicinity of the nucleus the 
brightest component is the one originating in the galactic 
disc, while away from the nucleus the brightest component 
is the outflowing one. As a consequence, the centroid veloc- 
ities show an apparent increase in velocity outwards, that 
mimics acceleration along the NLR. In the channel maps 
we see high velocity gas at the nucleus, showing that the 
outflow does not leave the nucleus at zero velocity. We in- 
terpret this high velocity gas as an outflowing wind from 
the AGN or ambient gas more directly interacting with this 
wind. As this highest velocity gas moves away from the nu- 
cleus, it pushes and accelerates the gas from the disc. The 
small amount of high velocity gas thus transfers momentum 
to a large reservoir of low velocity gas from the disc, caus- 
ing the observed emission-line profiles and velocity structure 
along the bi- cone. This scenario is co nsistent with the one 
proposed bv lEverett fc Murray! (|2007f ) who, after verifying 
the failure of various wind models to reproduce the observed 
outflows in NGC 4151, suggested that the interaction of an 
outflowing nuclear wind with the ISM could explain the ob- 
served velocity profile. 
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Figure 4. Velocity dispersions obtained from the fits of the single (top) and double (bottom) components for [Sin] (left), and He I 
(right). The continuous line shows the orientation of the galaxy major axis and the dashed line shows the orientation of the bi-cone. 
Emission-line profiles and fittings typical of the regions marked with letters A, and E are shown on Fig. [T] for the [S III] emission line. 



5.2 Relation between the gas kinematics and the 
radio jet 

Some of the previous studies of the NLR of NGC 4151 have 
suggested that there is n o relation between the radio jet 
and the bi-conical outflow jDas et al. l l2005l ; ICrenshaw et al.l 
l200d ; lHutchings et al.ll 19991 ). The distinct orientations of the 
two outflows do support distinct origins for them: the radio 
jet (oriented at PA« 75° ) originates in the funnel of the ac- 
cretion disc, while the bi-conical outflow originates further 
out, probably still in the accretion disc. The disc should also 
be warped, in order to produce a wind with a distinct ori- 
entation (PA=60°) from that of the radio jet. The channel 
maps of Figs. [5] and [7] confirm that the flux distributions 
at high velocities show no relation with the radio jet. But 
at velocities close to systemic the flux distributions become 
oriented along the radio axis instead of along the bi-cone. 
In addition, a comparison between the channel maps for 
velocities close to systemic and the contours of the radio 
map reveals line flux enhancements at the locations of ra- 
dio knots. Our proposed interpretation is that the radio jet 
has been launched close to the plane of the galaxy, and does 
interact with the circumnuclear gas pushing it out approx- 
imately along the plane and introducing the flux enhance- 
ments mostly due to gas density enhancement which are 
observed at velocities close to systemic. 

The observed flux enhancements could also be due to 
extra ionization by shock-induced X-rays in the NLR clouds 
interacting with the jet. Recent Chandra X- ray observations 
of the NLR of NGC 4151 IWang et al.l (|2009T ) have shown ex- 
tended emission along the NLR with X-ray flux enhance- 



ments in association with radio knots. They show that the 
radio-jet parameters are only consistent with the X-ray flux 
beeing due to thermal emission associated with the radio 
jet. 



The presence of shock heating is supported by the re- 
sults of Paper I, namely the increase of the [Fe n] 1.257/[P n] 
line ratio to values larger than ~ 2 in association with the 
radio knots. Values larger than « 2 imply the presence of 
fast shocks to destroy dust grains and release the Fe. We 
also find, in Paper I, a high temperature ( 16 000K) in the 
flux-enhanced regions, as derived from the [Fe n] line ratios, 
in agreement with this scenario. In the channel maps, the 
relation of the gas emission with the radio jet seems also 
to be strongest for [Fe n] , as can be observed in Fig. [6] at 
velocities from ~ —200 to ~ +100 kms -1 , there is enhanced 
[Fell] emission at the locations of most radio knots. There 
seems also to be an association between a "flared structure" 
to the SW observed in the radio map with the line emission 
at velocities of up to ~ —300 kms" 1 . This suggests that the 
radio jet has been deflected towards the line-of-sight at this 
location. Some spatial association between optical and radio 
emission in the NLR of NGC 4151 has also been observed by 
iMundell et all |2003). Comparing high resolution radio im- 
ages to an HST [O ill] image, they found a number of bright 
[O ill] clouds that seem to bound a number of knots observed 
in the radio image, which they interpret are due to shocks 
of the radio jet with the circumnuclear ISM as the jet moves 
through the NLR and encounters gas clouds. 
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Figure 5. Velocity channels obtained by integrating the flux within velocity bins of 63 kms~ 1 along the [S III] emission-line profile. The 
numbers in the upper left corner of each panel are the central velocity of the bin, in kms" 1 relative to systemic. The continuous line 
shows the orientation of the galaxy major axis and the dashed line shows the orientation of the bi-cone. Green contours are from the 
radio MERLIN image. 



5.3 Comparison with previous studies 



ICrenshaw fc Kraemerl (2000) studied the [O ill] kinematics 
of the NLR of the prototypical Seyfert 2 galaxy, NGC 1068, 
using a long-slit spectrum obtained at PA 202° with STIS 
on HST. They found a distribution of line centroid veloci- 
ties along the slit that resembled a "figure of eight", with 
prominent high- velocity redshifted and blueshifted emission 
clumps and relatively few low- velocity clouds at intermedi- 
ate distances from the nucleus. This distribution led them 
to infer a highly inclined conical geometry for the NLR of 
NGC 1068, with the emission arising from clouds outflow- 
ing radially along the cone walls. The cloud outflow veloci- 
ties increase approximately linearly with radial distance on 
both sides of the nucleus to ~ ±500 km s -1 at distances of 
~ ±2 arcsec fro m the nucleus, and then decli ne back to zero 
by ~ ±4 arcsec. ICrenshaw fc Kraemerl (2000) were forced to 
hypothesise a gradual acceleration of the NLR clouds over 
~ ±140pc followed by gradual dec eleration to ~ ±3 00 pc to 
account for this velocity structure. iDas et aL I (|2006l ) used a 
more extensive set of STIS spectra to develop this model for 
NGC 1068 further. 

ICrenshaw et al.l (|2000T ) applied the same model to STIS 
spectra of the NLR of NGC 4151. The NLR of NGC 4151 
is far less incline d to t he line of sight than NGC 1068 
- ICrenshaw et al.l (2000) infer 50° compared to 85° for 
NGC 1068 - so radial flows along the cone walls do not sepa- 
rate as cleanly into high-velocity redshifted and blueshifted 



structures. Consequently, the fit to the NGC 4151 d ata is less 
compelling than for NGC 1068. IDas et all (|2005h revisited 
the NGC 4151 NLR kinematics by fitting a similar model to 
spectra from five adjacent STIS slit positions. Their compre- 
hensive analysis provides strong evidence that the NLR of 
NGC 4151 is inclined by 45° to the line of sight, with a hol- 
low bi-conical outflow o c cupyi ng cone half-angles between 
15° and 33°. IDas et all (|2005T ) conclude that the outflow 
along the walls of the bi-cone is best reproduced with a ve- 
locity law: v oc r valid for distances r from the nucleus up 
to a turnover point at r = 96 pc (= l'.'48), from where the 
gas starts to decelerate and reaches the systemic velocity at 
400 pc (= 6'.' 15). The SW side of the bi-cone is tilted toward 
us, as illustrated in Fig-H] In this geometry and orientation, 
we view the nucleus just outside the outer wall, at an an- 
gle of 12° with respect to the nearest part of the SW cone. 
This geometry and orientation has led to the idea that we 
see NGC 4151 as a Seyfert 1 galaxy because the BLR can be 
viewed directly through a clumpy medium defining the cone 
walls. The angle betwe en the bi-cone ax is and the normal to 
the galaxy disk is 36° (|Das et al.l l2005). and thus, for a to- 
tal opening angle of the bi-cone of 66° , one side of the cone 
exits the galactic disk steeply, with the outer wall making 
an angle of 87° with the galaxy plane, while the other outer 
wall leaves the plane at only 21°, impacting a larger area 
of the disk (for a typical scale height for the disc of several 
tens of pc), as illustrated in Fig. [9] 
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Figure 7. Channel maps for the Br7 emission line, as in Fig. \E\ 
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Figure 8. Channel maps for the H2 emission line, as in Fig. [5] The dotted line shows the orientation of the galaxy minor axis. 
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Figure 9. Bicone geometry as proposed bv lDas et al ] |2005h . In 
the left panel we show a representation of the bi-cone and galaxy 
disk at the observed orientation on the sky, while in the right 
panel the bi-cone is shown edge-on. 



5.3.1 Centroid velocity maps 

IDas et all (|2005h argue that the intersection of the cone with 
the disc produces the observed geometry of the ex tended 



narrow-line region as shown by Evans et al.l l| 19951 ). If we 
now inspect our centroid velocity ma ps of Fig. [2 ] and F ig. [3] 
in light of the geometry proposed by IDas et al. |2005l ). we 
conclude that, in most regions where we see two velocity 
components, there is always one component at the systemic 
velocity plus an outflowing component in blueshift to the 
SW and in redshift to the NE. We attribute the systemic 
velocity component to emitting gas from the galaxy disc, 
in which case we are led to the conclusion that the AGN 
ionizes a region of the galaxy disc comparable in extent to 
that of the outflowing gas emission we detect along the bi- 
cone axis (~l'.'5, or «100pc) and somewhat less extended 
perpendicular to it (~1", or ~65pc). 

Blueshifted components are observed to the SW at ~ 
— 250 kms -1 and at higher velocities of up to —600 kms -1 
(Fig. [2] for [S in] and Fig. [3}. We interpret the component 
with blueshifts ~ —250 kms -1 at the position labelled D in 
Fig. [2] together with the single component surrounding it 
(region iJin Fig. [2}, as originating in the back surface of the 



SW cone. In the bottom part of the region we labelled C in 
Fig. [2] we observe a few points with higher blueshifts of up 
to —600 kms -1 (see also Fig|3J. These components probably 
originate in the front part of the SW cone, which we see at a 
small angle (12°) and thus with a large velocity component 
along the line-of-sight. There is less emission from this front 
part of the SW cone than from the back, and our hypothesis 
to explain this is that the cone walls are formed by gas 
entrained from the galaxy disk. The front wall is thus likely 
to be comprised of more tenuous gas than the back wall, as 
there is less gas to be entrained at high angles to the galactic 
disk. 

Redshifted components are observed to the NE at some- 
what lower velocities than to the SW (region labelled B in 
Fig. [2}. This can be seen more clearly in Fig. [3] there is a 
sequence of redshifts of up to ~ 150 kms -1 , which can be 
attributed to the front of the NE cone and then some higher 
redshifts of up to 400-500 kms -1 , which can be attributed 
to the back of the NE cone. In Fig. [2] there is more emission 
from gas at the highest redshifts than emission from gas at 
the highest blueshifts. This result, together with the lower 
redshifts to the NE than blueshifts to the SW, suggest that 
the NE cone has its axis somewhat tilted towards the line 
of sight, instead of being aligned with the axis of the SW 
cone. This tilt would allow us to see more of the back wall 
of the NE cone and would result in lower observed velocities 
for both its front and back walls. 

In the model proposed by IDas et ail (120051 ). the NLR 
gas first accelerates up to ~ 100 pc from the nucleus and 
then decelerates, analogous to their successful model for 
NGC 1068. While it is appealing in terms of AGN unifica- 
tion for NGC 1068 and NGC 4151 to be fit with a similar 
model, neither our centroid velocity maps of Figs. [2] and [3] 
nor the [O ill] kinematic data for NGC 4151 unambiguously 
point to this interpretation. The gradual acceleration over 
the inner ~ 100 pc and th en deceleration required by the 
model remain unexplained. lEverett fc Murray! i|2007l ) noted 
that radiative and magnetic acceleration from the region of 
the accretion disk both accelerate gas on length scales of the 
order of the launch radius, so winds driven by these mech- 
anisms should reach their terminal velocities by ~ 10 pc 
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from the central source, not the ~ 100 pc required by the 
mode l fits to the data. Therm al winds also fail in this re- 
gard ijEverett fe Murravl l2007'l . Dynamical models including 
various acc eleration and deceleration mechani sms have been 
developed l|Das. Crenshaw, fe Kraemer] I2007T ) . While these 
models fail to identify the acceleration mechanism, they do 
demonstrate that gravitational forces alone cannot produce 
the required deceleration, while drag forces from the inter- 
action of NLR clouds with a diffuse, hot, X-ray-emitting 
ambient medium could be responsible. 

The difficulty in identifying a suitable wind accelera- 
tion mechanism leads to the question of whether some of 
the NLR kinematics are due to gradual global wind accel- 
eration or whether they result from the local interaction of 
an already accelerated wind with its ambient medium, as we 
have proposed above on the basi s of the channel ma p s. Thi s 
scenario has been proposed by lEverett fc Murravl (|2007h . 
who realised that the interaction of an outflowing nuclear 
wind with the ISM could reproduce the apparent accelera- 
tion o f the NLR following the works of iMatzner fc McKeel 
l|l999l ) and lLada fc Fichl |l996l ). These authors showed that 
a protostellar hydromagnetic wind expanding into the sur- 
rounding gas produces a bipolar outflow with velocity pro- 
portional to the distance to the star until a maximum radius 
beyond which the gas starts to decelerate. 



5.3.2 Position-Velocity diagrams 

In order to be abl e to compare our da t a with thos e of pr e- 
vious works (e.g. iKaiser et all (|2000h . iDas et all (|2005h ). 
we have built position-velocity diagrams (hereafter PV dia- 
grams), which are shown in Fig. [10] for the [S ill] emission 
line. These diagrams were obtained from the Z-band NIFS 
data cube by extracting 0.2 arcsec wide pseudo-slit spectra 
oriented at PA= 60° (which is the orientation of the bi-cone 
axis), offset from the nucleus perpendicular to the pseudo- 
slit direction by between -0'.'8 (along PA= 150°) and +0'.'8 
(along PA= -30°) in steps of 0'.'2. The [S in] diagrams are 
plotted both on a logarithmic scale (top of Fig. 1 10 p to reveal 
the faint emission structure and on a linear scale (bottom) 
to show the regions with strongest emission. While qualita- 
tivel y similar to pre vious long-slit [O ill] spectra (e.g. Fig. 
2 of Das et al"1l2005h . the [Sin] diagrams seem to show less 
emission at velocities exceeding ±1000 k ms" 1 . In par t icular , 
high-velocity [Om] clouds A and B in IKaiser et al.l (120001 ) 
(near — 1400 kms -1 and — l'.'l) are not apparent in our data. 
In part, this could be due to the poorer spatial resolution 
obtained with NIFS+ALTAIR at short near-infrared wave- 
lengths. 

As with the [O ill] line, most of the [S ill] emission occurs 
at velocities below ~ ± 700 kms . The central diagrams of 
Fig. [10] can be compar ed directly with the [O ill] slit spec- 
trum shown in Fig. 2 of IDas et al.l (2005). The [S ill] spectra 
are much less structured than the [O ill] spectra, presum- 
ably due to poorer spatial resolution and consequent lower 
sensitivity to faint emission. Nevertheless, the two emission 
lines show a similar overall appearance along the central slit 
with blueshifted emission near —250 kms -1 extending from 

~ -3" to r 0'.'3 from the nucleus towards PA= -120° 

(SW, bottom of Fig. I10p . emission at the systemic veloc- 
ity near the nucleus, and some redshifted emission near 
+ 150 — 200 kms" 1 extending to ~ 1 arcsec from the nu- 



cleus towards PA= 60° (NE, top of Fig.[10]). The off-nucleus 
PV diagrams in Fig. 1 101 have a similar overall structure, but 
show that the — 250 km s -1 component is present in all neg- 
ative offsets up to the slit at — 0'.'8. Of note is the fact that 
the velocity of this emission remains approximately constant 
over its full extent. The brightest +200 kms -1 emission is 
only clearly apparent towards PA= 60° along the cone axis 
(i.e., in the upper part of the central panel of Fig. I10[) . 
Some aditional emission from this component can be ob- 
served along the slits with positive offsets at fainter levels 
and at larger distances from the nucleus. The systemic ve- 
locity component is prominent in the upper section (NE) of 
pseudo-slits with negative offsets and in the lower sections 
(SW) of pseudo-slits having positive offsets, and originates 
in circumnuclear gas from the galaxy plane. 

Fig. [TTJ shows the corresponding PV diagrams for the 
[Fell] 1.644 (i.m line. The emission in the central pseudo- 
slit corr esponds closely t o the brightest [O ill] emission in 
Fig. 2 of IDas et al.l (|2005T ). but the high- velocity extensions 
of the [O ill] clumps are much less prominent in our [Fe n] 
data. Although the major velocity components apparent in 
the [S ill] line are also apparent in [Fe n] , the higher spa- 
tial resolution achieved at this longer wavelength helps to 
reveal sharper structures than in Fig. 1101 On the other 
hand, the structures seen in the [Fell] emission can also 
in part be attributed to enhancements associated with the 
radio jet, as observed in the channel maps. The compo- 
nent near —250 kms -1 is comprised of several clumps and 
extends to only ~ — 1'/5 along the central pseudo-slit to- 
wards PA= —120°, but reaches —2 arcsec at large nega- 
tive, and possibly also positive, slit offsets - i.e., in the ex- 
treme left and right panels of Fig. 1111 Clouds with velocities 
« 250 kms -1 are observed at ~ +1 arcsec along the central 
slit and closer to the nucleus along the 0'.'2 and 0'.'4 offset 
slits. Emission with velocities up to +500 km s -1 is appar- 
ent at slit offsets of — 0'.'2 and +0'.'2. The systemic-velocity 
component is restricted to the nucleus in the central slit, but 
extends to nearly —1 arcsec along the +C/.'2 to +0'.'6 offset 
slits and to nearly +1 arcsec along the — 0'.'2 and — 0'.'4 off- 
set slits. However, the systemic- velocity emission at negative 
slit offsets - to the left in Fig. [TTJ- is dominated by a further 
bright clump near +1 arcsec - approximately east of nucleus, 
in the vicinity of the receding radio lobe. Overall the impres- 
sion from the [Fe n] data is of roughly constant redshifted ve- 
locities towards PA= 60° and blueshifted velocities towards 
PA= —120° with an additional systemic- velocity component 
originating in the galaxy disc and associated with the radio 
jet. 



5.4 Kinematic models for the NLR 

The role of the geometry and projection effects on the obser- 
vations can best be evaluated via modelling. In this section, 
we compare our data with models with and without accel- 
eration along the NLR. 



5.4-1 Das et al. Model 

We begin by testing the model proposed bv lDas et al.l (|2005T l 
on the basis of the optical [O ill] emission obtained with 
STIS, namely of an inclined hollow cone with radial out- 
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Figure 10. Position-velocity diagrams of [S III] 0.9533 um emission for 0.2" wide pseudo-slits oriented at PA=60° along the cone axis 
and offset perpendicular to the pseudo-slit by -0.8" to +0.8" as indicated in each frame. An average continuum image extracted between 
±(2000-2500) kms - 1 has been subtracted. Crosses mark the location of the nucleus in the central pseudo-slit and the "equivalent" 
location for the other slits. Top: Line intensities are displayed on a logarithmic scale between 4.6 X 10" 
A -1 arcsec" 2 . Bottom: Line intensities are displayed on a linear scale. 
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Figure 11. Position-velocity diagrams of [Fell] 1.6440 |om emission shown as in Fig llOl 



flow along the cone walls that is characteri sed by linear ac- 
celera tion followed by linear deceleration (|Crenshaw et al.l 
2000). Our 3D kinematic data allow us to test this model 
against both channel maps and PV diagrams through multi- 
ple pseudo-slits in ways that were not available to previous 
authors. We begin b y presenting bot h these views for the 
model parameters of iDas et alj (|2005r ) in Figs. 1121 and [131 



As our data cover essentially the entire region of accelera- 
tion in the models, emission from decelerating gas beyond 
the velocity maximum has been suppressed in Fig 1131 



The IDas et ail (|2005l ) model predictions differ in sev- 
eral ways from the near-infrared emission-line data. First, 
the high velocity gas just to the SW and NE of the nucleus 
predicted by the models (Fig. I12[) and seen in [O ill] is less 
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Figure 12. Simulated position- velocity diagrams for the hollow cone with radial outflow model of lDas et al.l l|2005h . The orientation and 
pseudo-slit geometry are as for Figs. ITOl and 1 11 1 
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Figure 13. Simulated channel maps for the hollow cone with 
radial outflow model of iDas et all (|2005h . The orientation and 
velocity bins are as for Fig.0 



conspicuous in our data. It is seen in our data only in the 
log plots for [S in] , where the faint emission gets enhanced 
(top panel of Fig. [10} . We also see little evidence of accelera- 
tion and deceleration in the near-infrared data, which mostly 
show sequences at constant velocities from the nucleus all 
the way out. Even in the deep [S in] PV diagrams, gas up to 
3 arcsec from the nucleus in the lower half of these diagrams 
- extending to PA= —120° - is still moving at blueshifted 
velocities of ~ — 250 km s -1 , with similar but smaller red- 
shifted velocities in the upper half of the diagrams. 

Another significant difference between the model and 
the near-infrared PV data is that the model predicts that 
pseudo-slits offset from the nucleus will show no emission 
near the slit centre - Fig. 1121 This is the equatorial region 
away from the nucleus and between the two oppositely- 
directed cone "lobes". Most of the actual emission at off- 
nucleus slit positions arises in this region that is excluded 
by the model - c.f., Fig. 1111 There is also emission at zero 
velocity - mainly observed in the [S ill] PV diagram which 
is not included in the model. 

Finally, the model predicts that the channel maps 
should show ridges of emission perpendicular to the cone 
axis at high positive and negative velocities - Fig. 1131 This 
morphology arises because outflowing material in the hollow 
cone essentially emits from an annular region perpendicular 
to the cone axis - convolved with other orientation factors. 



In contrast, the channel maps in Figs. [5] and [6] show that 
the highest-velocity emission is confined to positions close 
to the cone axis, and is possibly extended along the cone 
axis towards the nucleus. 



5.4.2 Modified Das et al. Model 

We now investigate modifications to the hollow-cone, radial- 
outflow model that attempt to address the above differences. 
The simplest modification is to hypothesise that the cone 
wall that is steeply inclined to the plane of the sky does not 
emit strongly in near-infrared lines. This could occur be- 
cause this cone wall is more inclined to the galactic disk. As 
discussed in Sec. 15.3.11 under the hypothesis that the cone 
walls are formed by gas entrained from the galaxy disk, this 
wall is likely to be comprised of more tenuous gas, as there 
is less gas to be entrained at high angles to the galactic disk. 
This modification removes the high-velocity components in 
Fig. 1121 To this, we then add an illustrative systemic-velocity 
component that is aligned with the radio jet and an arbi- 
trary galactic disk component. The radio-jet component is 
modeled as a central nuclear component plus two clumps 
offset E-W from the nucleus by 1". Each clump is assumed 
to be at the systemic velocity of the galaxy. The galaxy com- 
ponent has a Gaussian intensity profile and is assumed to be 
at the systemic velocity. The presence of this component is 
suggested both by the zero velocity component in the cen- 
troid velocity maps - Fig. [2] - and by the channel maps, 
which show strong emission at zero velocities (Figs. [5] and 
[6}. The resulting PV diagrams are shown in Fig. 1141 

The PV diagrams in Fig. Q3] do not provide a perfect fit 
to the data (in particular to the highest spatial resolution 
[Fell] data in Fig. ll 1 1) . However, they serve to illustrate the 
following points: Firstly, the emission from the cone (blue 
regions in Fig, ll4[) is a better fit to the data in the linear 
plots (showing the strongest emitting gas) if only the sec- 
tion of the cone wall closest to the plane of the sky, and 
hence closest to the galactic disk, is included. This may be 
due to interaction of the NLR gas in the cone wall with the 
galactic disk, possibly by the entrainment of disk material 
into the NLR outflow. Given that this material is moving in 
a direction closer to the plane of the sky, we are less sensitive 
to its true velocity so cannot make strong claims about the 
length scale on which it is accelerated. Fig. [14] is drawn for 
the gradual acceleration and subsequent deceleration pro- 
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Figure 14. Simulated position-velocity diagrams for the hollow cone with radial outflow model of lDas et al . (2005) modified to include 
only the hemisphere of the cone with lowest inclination to the plane of the sky. The outflowing gas emission is shown in blue, an E-W radio 
jet component is shown in red, and an arbitrary galactic disk component is shown in green. The orientation and pseudo-slit geometry 
are as for Fig. 1101 
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Figure 15. Simulated channel maps f or the hollow cone with 
radial outflow model of lDas et al.l (120051 ) modified to include only 
the side of the cone with lowest inclination to the plane of the 
sky. The orientation and velocity bins are as for Fig. 1101 

posed by iDas et ail (|2005h . However, the data permit vir- 
tually any velocity law, including constant velocity outflow, 
on the length scales to which they are sensitive. The second 
point illustrated by Fig. [14] is that emission from gas inter- 
action with the radio jet is significant. This is modeled as 
nuclear and off-nuclear clumps, but in reality other emit- 
ting gas forms a bridge between the main clumps and the 
nucleus. This is seen in the channel maps shown in Fig.|51 Fi- 
nally, the galaxy component at systemic velocity reproduces 
the data at PA=60° for pseudo-slits at negative offsets and 
at PA= —120° for pseudo-slits at positive offsets, supporting 
the presence of a disc component, as concluded also from 
the centroid velocity maps of Fig. [2] 

Simulated channel maps are shown in Fig. 1151 These 
only approximately reproduce the data presented in Fig. [6] 
In particular, two separated emission clumps aligned per- 
pendicular to the cone axis are predicted from the cone walls 
at high-redshifted and high-blueshifted velocities while only 
a single clump near the cone axis is actually seen in the data. 

5.4-3 Constant Velocity Model 

Finally, we have considered a truncated hollow conical out- 
flow at a constant velocity of 600 km s -1 , as we could not 



see acceleration or deceleration in our data. As in the pre- 
vious model, we add the radio jet and galactic disk compo- 
nent. The resulting position- velocity diagrams are shown in 

Fig. [E! 

The fit to the observed constant velocity sequence at 
~ — 250kms _1 for [S ill] (Fig. [TO]) along the slit passing 
through the nucleus and the ones with negative offsets are 
improved by this model. The disc component at zero veloc- 
ity reproduces emission observed at the nucleus and along 
PA= 60 deg (positive position in the panels of Fig. 1 10 p and 
PA= 150 deg (to the left in Fig. 1 1 p , as well as less extended 
emission along PA= 240 deg (negative positions in Fig. 1 1 p 
and PA= 330 deg (to the right in Fig.flfJ). The high veloc- 
ity sequences predicted by the models - although not con- 
spicuous in the linear PV diagrams of [S in] (bottom panel 
of Fig. I10p - are observable in the log PV diagrams (up- 
per panel of Fig. llOp . due to its faint emission. These high 
velocity sequences and are due to the front wall of the ap- 
proaching cone and far wall of the receeding cone. We point 
out that this short high velocity sequence as predicted by 
the constant velocity model seems to be adequate to repro- 
duce also t he hig h velocity gas observed in [O in] (Fig. 2 of 
IDas et al.l (|2005l )V 

Nevertheless, this model did not improve the fit for 
pseudo-slits offset by more than (f.'A from the nucleus, which 
still do not show emission near the slit centre in the model. 
In addition, these pseudo-slits show velocity sequences that 
decrease with increasing distance from the nucleus. In or- 
der to improve the fit, new geometries should be considered. 
For example, a flow that is not conical, but in the form of a 
bowl as observed in planet ary nebulae and in re cent NIFS 
observations of NGC 1068 (jMcGregor et alj|2010h . 

Simulated channel maps are shown in Fig. 1171 
These maps are an improvement relative to the previ- 
ous model for velocity channels between — 300kms _1 and 
+300kms _1 , but show double structures for velocities be- 
tween 300kms _1 and 500 kms , which are not seen in the 
data. For the highest velocities, a small triangular structure 
is seen, while the data show only clumps near the cone axis. 
This discrepancy can also be attributed to possible devia- 
tions from a bi-conical geometry and to very weak emission 
from the wall of the bi-cone that is most inclined. 
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Figure 16. Simulated position-velocity diagrams for the truncated hollow cone model with constant velocity 600 km s -1 outflow. The 
outflowing cone emission is shown in blue, an E-W radio jet component is shown in red, and the galactic disc component is shown in 
green (see text). The orientation and pseudo-slit geometry are as for Figs. [TOland llll 
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Figure 19. Integrated [Fell] 1.644 urn intensity formed by sum- 
ming the velocity slices in Fig. \E\ 



Figure 17. Simulated channel maps for the truncated hollow 
cone model with radial outflow at constant velocity, as in Fig. 1141 
(see text). The orientation and velocity bins are as for Fig. 1101 
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Figure 18. Simulated integr ated intensi t y ma p for the modified 
hollow cone outflow model of iDas et al.l l|2005h (see text), to be 
compared to Fig. 1191 



5.4-4 Integrated Intensity 

The modelled integrated emission-line map is similar for the 
three models discussed ab ove. We show in Fig. [18] the map 
obtained for the modified IDas et all (|2005l ) model. A com- 
parison with Fig. 1191 shows that it provides a reasonable re- 
production of the data, taking into account that one should 
add the contribution of gas from the galaxy plane to the 
modelled emission. 



5-4-5 The origin of the highest velocity gas seen in [Olll] 

ICecil et all (|2002h identified blueshifted [O in] emission with 
velocities up to 3200kms _1 in the NGC 1068 NLR 70 - 
150 pc NE of the nucleus. This is the region attributed to 
the highest radial outflow velocities in the gradual acceler- 
ation model. Many of these knots are linked kinematically 
to more massi ve and slower-mo ving clumps located closer 
to the nucleus. ICecil et al.l (|2002h identify these knots as ab- 
lation flows from disintegrating molecular clouds that are 
being photoionis e d and radiatively accelerated by the AGN. 
iHutchings et all l|l999l ) detected faint, high- velocity [O ill] 
emission with radial velocities up to 1400kms _1 in STIS 
slitless spectra of NGC 4151. The high-velocity emission was 
shown to arise in high-excitation ionise d gas, and has veloc - 
ity dispersions in excess of 130 kms -1 (|Kaiser et aLl feoOO), 
which is also consistent with an origin in an ablation flow. 

Near-infrared [S ill] emission arises in conditions similar 
to the optical [O ill] emission, and can thus be formed in 
tenuous fully ionised gas experiencing a high local ionisation 
parameter. But then, why do we see less [S ill] than [O ill] 
emission? Either we lack the sensitivity to detect faint high- 
velocity emission, or this emission is not present. 

In order to investigate this, versio n 8.0 of the photoion - 
isation code Cloudy, last described bv lFerland et al.l (|l998l ). 
was used to calculate the ratio of [O in] to [S in] for NLR 
clouds with a range of densities. The ionisation parameter 
was set to log U = —2.0 for a gas density of 500 cm -3 and 
then varied inversely with the gas density to simulate clouds 
with different densities being illuminated by the same radi- 
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log U 

Figure 20. [O III] / [S III] line ratios for different ionisation param- 
eters or densities with a fixed radiation field that corresponds to 
an ionisation parameter of log U = —2.0 at a NLR cloud with a 
density of 500 cm -3 . 



ation field. Such a situation would occur if the low-velocity 
[O ill] is emitted at the front face of a dense NLR cloud and 
the high- velocity [O ill] emission arises in an associated low- 
density ablation flow. The resulting emission-line ratios are 
plotted in Fig. [201 where it can be seen that [O m]/[S ill] in- 
creases from ~ 3 at densities typical of NLR clouds to ~ 13 
at densities < 10 cm -3 that might be typical of an ablation 
flow. The values at high densities are typical of the inte- 
grated line ratios observed in Seyfert galax ies NLRs (e.g., 
IStorchi-Bergmann. Kinney, fc Challisl [1995) . The strength 
of [S in] falls sharply to lower densities while the strength 
of [O ill] is largely unchanged. The reason for this depen- 
dence on ionisation parameter is that sulphur exists pre- 
dominantly as [S in] at typical NLR densities, but as [S iv] 
and higher species in more highly excited low-density re- 
gions. In contrast, the higher ionisation potential of [Oiii] 
(54.9 eV compared to 35.1 eV for [S hi] and 47.2 eV for 
[Siv]) means that [O in] survives better in highly excited 
regions. The lower intensity of [S ill] in low-density irradi- 
ated regions makes it more difficult to detect a low-density 
ablation flow in [S ill] emission than in [O ill] emission. This 
combined with the lower Strehl ratio and consequent poorer 
image quality achieved by adaptive optics systems at short 
near-infrared wavelengths probably account for our inability 
to detect more high-velocity gas in NGC4151. 

There is little evidence in our near-infrared data for 
the gradual acceleration follo wed by gradual deceleration 
proposed by iDas et all (|2005l ) based on [O in] kinematics. 
No viable physical mechanism has been identified for such 
a gradual acceleration. A possible explanation is that the 
high- velocity [O ill] emission does not represent the bulk 
motion of the outflow. Rather, this could be low-density 
material ablated from NLR clouds that we observe outflow- 
ing at »; 250 kms and are accelerated much closer to the 
AGN. This is consistent with the low surface brightness of 
the high-velocity, low-density emission in [S ill] and with the 
prominence of near-infrared emission formed on the sides of 
the outflow closer to high-density gas in the galactic disk. 



While the kinematics of the [S ill] , He I, Hi, and [Fell] emit- 
ting gas is dominated by outflows along the bi-cone, the 
H2 radial velocity field in Fig. [8] is dominated by velocities 
close to zero around the galaxy minor axis. Although the 
velocity range is small, it can be observed that, from the 
blueshifted to the redshifted channels, the flux distribution 
moves around the minor axis of the galaxy, from the SW to 
the NE, as expected for gas rotating in the galaxy plane. Al- 
though we see a kinematic component at systemic velocity 
also in the ionised gas, we do not see the outflowing compo- 
nent in H2. The velocity dispersion in the H2 gas is also the 
smallest among the emission-line profiles. 

Distinct kinematics for the H2 emitting gas relative to 
the emission lines of ionised gas h as also been observed in 
previous studies of other objects (IStorchi-Bergmann et all 
Il999l ; iRodrfguez-Ardila et alj l2005l : iRiffel et alJ|20Q6j ). The 
H2 emission lines usually have small velocity dispersions and 
their kinematics seem to be dominated by rotation in the 
galaxy plane, while the other emission lines have larger ve- 
locity dispersions and other velocity components besides ro- 
tation, as observed also in NGC4151. In Paper I we have 
argued that the flux distributions, physical properties and 
location suggest that the H2 emission gas traces the AGN 
feeding, while the ionised gas traces its feedback. Although 
we do not see the H2 emitting gas actually flowing down into 
the nucleus, it does present a distinct kinematic signature 
relative to the ionised gas. The inflow may be occurring fur- 
ther in, in non-emitting colder molecular gas. The hot molec- 
ular gas is estimated to be a fract ion of only 10~ r to 10~ 5 of 
the total mass of molecular gas (|Dale et al.ll2005h in AGN. 
In the case of NGC 4151, for which we obtained a mass of 
hot H2 of w 240 Mq (|Paper J ), we can thus estimate a total 
molecular gas mass of ~ 10 8 Mq. Such a large molecular gas 
reservoir around the nucleus of active galaxies is supported 
by radio observations of the central region of active galaxies 
in general, by, for example, the "NUGA" group, which re- 
port molecular gas masses in the range 10 7 -10 9 M B in the m- 



ner few hundred parsecs of active galaxies (e. g. Boone et all 



120071 ; iGarcia-Burillo et alj|2005 iKrips et al.N2007h . If such 
a reservoir of molecular gas has accumulated by a flow that 
was the trigger of the present nuclear activity, it can be 
concluded that a steady flow would have to bring gas in- 
wards at a rate in the range 0.1-10 M Q yr _1 . In NGC 4151, 
iMundell fc Shonel l| 19991 ) have measured inflows in radio ob- 
servations of Hi that occur along narrow channels inter- 
preted as shocks along the leading edges of the large scale 
bar (at PA= 130°). These inflows seem to end at the inner 
part of the bar in the nuclear region. The orientation of the 
flow is approximately perpendicular to the "arcs" seen in the 
H2 flux distribution - see Fig. [5]- and one possibility is that 
the radio observations are tracing the gas inflows that accu- 
mulate in the nuclear region and contribute to build up the 
molecular gas, part of which we observe due to excitation of 
the H2 molecule by X-rays emitted by the nucleus. 

The large amounts of molecular gas observed around 
the nuclei of active galaxies can only lead to the conclu- 
sion that inflows must be occurring |Garcia-Burillo et all 
2005). Nevertheless, the actual kinematic observation of in- 
flowing gas in active galaxies is difficult to obtain. This is 
due to the fact that the emission of the cold, inflowing gas 
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is much w eaker than the em issio n of the hotter, outflow- 
ing ga s. In lFathi et al.l l|2006l ) and IStorchi-Bergmann et alj 
(2007) we were successful in mapping the inflow in low- 
excitation ionised gas because the targets were low-activity 
AGN - the LINERs NGC 1097 and NGC 6951 - for which 
the nuclear radiation and outflows are much weaker than 
in Seyfert galaxies, thus allowing the observation of inflows 
in the disc in ionised gas. We are probably seeing just the 
ionised "skin" of a flow of cold gas on its way to feed the 
AGN. We were also successful in mapping an infl ow observed 
in ho t molecular gas in the galaxy NGC 4051 (I Riff el et al.l 
2008). In this latter case, the X-ray emission of the AGN 
may have worked as a flashlight which "illuminated" the in- 
flowing molecular gas. We obtained a mass inflow rate of 
4 x 10 -5 Mq yr _1 . If we use the same ratio of hot to to- 
tal molecular gas as above, the inflow rate becomes of the 
order of a few M^yr^ 1 . Another recent case of i nflow seen 
in molecular gas was reported bv I Sanchez et all (l2009h . us- 
ing near-IR IFU observations of NGC 1068 obtained with 
SINFONI. 

The feedback from the AGN - in the form of the mass 
outflow rate along the NLR and its kinetic power - can be 
calculated using the ionised gas luminosity and geometry of 
the outflow through the expression: 



M H n = m p Ami v f A 



(1) 



where m p is the mass of the proton, Ami is the ionised gas 
density, v is the outflow velocity, / is the gas filling factor 
and A is the cross section area of the outflow. The filling 
factor was estimated from: 



Lbtj = jBr-iV f 



(2) 



where Lsr-y is the Br 7 luminosity emitted by a volume 



V and jsr 



6.94 x 10 23 erg cm 3 s 1 for a tempera- 



500 cul- 



ture T=10000K and gas density Ami 
|Osterbroklll989l ).FI 

The vol ume V was calcu lated assuming the geometry 
proposed by iDas et all (|2005l ) and used in the modelling 
described in the previous section, namely a hollow bi-cone 
with a height of 130 pc for each cone (corresponding to 2"), 
and inner and outer opening angles of 15 deg and 33 deg, 
respectively. We obtained V = 2.88 x 10 61 cm~ 3 . The Br 7 
luminosity was obtained from Paper I. We have further di- 
vided this by a factor of 2, because part of the ionised gas 
emission comes from galactic disc gas that is not outflowing 
(see previous section). We estimate from the channel maps 
that approximately half of the total ionised gas luminosity 
arises in this galactic disc gas. Substituting these estimates 
in the expression above we obtain a filling factor value of 
/ = 0.11. Filling factor values reported in the literature 
are usually smaller, but they usually consider an isotropic 
(spherical) outflow. Assuming an isotropic outflow instead 
of a bi-conical one, we obtain / = 0.025 which is similar to 
the typical values in the literature. 

The cross section area A of the hollow cone was cal- 
culated at a distance of 65 pc (1") from the cone apex for 
the hollow cone geometry and resulted in A = 5 x 10 40 cm 2 . 



2 We have tried a lower density of 100 cm -3 as in |Paper"ll but got 
as a result f>l, which is unphysical, supporting a higher average 
gas density. 



Adopting the velocity v = 600 kms" 1 for the outflow (which 
is the one assumed for the constant velocity outflow model) , 
Ami = 500cm -3 , / = 0.11 we finally obtain a mass outflow 
rate of Mhii ~ 1.2Moyr -1 in each cone. This value is of 
the order of the mass outflow rates reporte d for the NLR 
of ac t ive galaxies in n eutral and ionised gas ()Veilleux et alJ 
120051 : iHolt et alj [20061 ). a nd ^6 times larger than the mass 
outflow rate estimated bv lCrenshaw fc Kraemerl l|2007l ) from 
UV absorbers in the spectra. 

The mass outflow rate can be compared with the accre- 
tion rate to feed the supermassive black hole. The accretion 
rate can be estimated from: 



Ltoi 



(3) 



where c is the speed of lig ht, and the bolomet ric luminosity is 
L bol = 7.3 x 10 43 ergs" 1 (|Kaspi et al.ll2005l ). which is 0.012 
times the Ed dington luminosity for a black hole mass of 



4.5 x 10' M l|Onken et a l. 2007; Bentz et al.ll2006l ). For this 



Eddington ratio we adopt an efficiency factor for the conver- 
sion of rest mass energy into radiation of 77 = 0.1, a typical 
value for an optica lly thin geometrically thick accretion disc 
jFrank et al.ll2002l ). The result is m = 1.3 x 10~ 2 M yr _1 , 
which is two orders of magnitude smaller than the mass out- 
flow rate, indicating that the origin of the outflowing gas is 
not the AGN but could be surrouding gas from the galaxy 
ISM accelerated by a nuclear outflow, as discussed in the 
previous section. 

Our data also allows us to estimate the kinetic power of 
the outflow, including both the radial and turbulent compo- 
nents of the gas motion. Assuming that the velocity disper- 
sion of the Hn emitting gas is due to turbulence, the kinetic 
power can be estimated as: 



E ! 



M, 



(v 



2 . 2n 



(4) 



From Fig. [4] we obtain a ~ 150kms _1 , and adopting as 
above, v = 600kms~ 1 and using Mhii ~ 2.4 Mg yr" 1 
(considering the outflows in both cones), we obtain E ~ 
2.4 x 10 41 ergs~ 1 , which is « 0.03xU oi . Thus, from the 
power that can be extracted from the accretion rate of 
m = 1.3 x 10 -2 A/q yr -1 , ~0.3% is transformed in kinetic 
power in the NLR outflows. We note that this fraction is 
of the order of the AGN feedback derived in simulations 
for co-evolution of black ho les and galaxies (0.5%) , in order 
to match the M-ct relation l|Di Matteo et alj[2005r ). This ki- 
netic power is somewhat larger than those we have recently 
derived for a s ample of 6 nearby S eyfert galaxies showing 
NLR outflows (|Barbosa et al.ll2009l ). which are in the range 
10" 5 < E/L bol < 10" 4 . 

In order to gauge the impact that the above kinetic 
power may have in the galaxy, we estimate the total ki- 
netic energy released if the outflow as now observed lasts 
for at least 10 r yrs (as an estimated duration for the ac- 
tivity phase of the galaxy) to be 7.6 x 10 55 erg. This en- 
ergy can then be equated to the binding energy of the gas 
mass that could be pushed away from the centra l regio n 
of the galaxy (following the reasoning in lHolt et~ al. (2006)) 
E bind « GM gas M tota i/R gaa . For E bind = 7.6 x 10 55 erg, 
Mtotai = x 10 9 Mq (estimated from the galaxy absolute 
magnitude) and R ga s ~ 500 pc (just beyond the observed 
outflow), we obtain M gas « 10 8 Mq, which is ~ 5 times 
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the mass which outflows through the NLR in 10 7 yrs. It is 
interesting to note that this mass is of the order of the es- 
timated molecular gas mass present in the nuclear region of 
NGC4151, as discussed above, being also the typical value 
observed in the nuclear region of active galax ies in general 
l|Garcfa-Burillo et al.ll2005l : iBoone et alJl2007t >. Thus, it ap- 
pears that in the case of NGC 4151 the molecular gas reser- 
voir that is needed to trigger the nuclear activity can be 
blown away by the NLR outflow by the end of the activity 
cycle. 



7 CONCLUSIONS 

We have measured the kinematics of the ionised and molecu- 
lar gas surrounding the nucleus of NGC 4151 using integral- 
field spectroscopy obtained with the NIFS instrument at the 
Gemini North telescope. The main results of this paper are: 

• The ionised-gas kinematics are dominated by three ve- 
locity components: (1) extended emission at systemic veloc- 
ity observed in a circular region around the nucleus; (2) a 
component outflowing along the bi-cone (PA= 60deg); (3) 
another component due to the interaction of the radio jet 
with the galactic disk. 

• Regarding the extended emission at systemic velocity: 
The origin of this component is gas from the galaxy plane 
ionised by the AGN up to at least 1 arcsec (65 pc) from the 
nucleus, revealing that there is escape of ionizing radiation 
even along the equatorial plane of the bi-cone. 

• Regarding the outflowing component: Channel maps 
extracted along the emission-line profiles show high- velocity 
gas all the way from the nucleus to the border of the emit- 
ting region. This seems to contradict previous studies show- 
ing acceleration along the NLR. Our data suggest that the 
NLR clouds are accelerated very close to the nucleus (within 
« 10 pc), after which the flow moves at essentially constant 
velocity. 

• We do not detect NIR counterparts to the highest- 
velocity [O ill] emission seen by previous authors. It may 
originate in very low-density gas, which could be ablated 
material flowing off NLR clouds. This emission would not 
be characteristic of the bulk of the NLR mass. 

• The origin of the bulk of NLR mass seems to be en- 
trained gas from the galaxy plane. This interpretation is 
supported by the large mass-outflow rate in ionised gas of 
~ 1 Mq yr _1 obtained for the outflowing component in each 
cone, which is ~100 times the nuclear mass accretion rate. 

• The estimated kinetic power of the outflow is « 2.4 x 
10 41 ergs _1 , which is ~0.3% of L(, (, and is » 10 times 
larger than values obtained previously by our group for other 
nearby Seyfert galaxies. 

• Regarding the third kinematic component: Our data 
clearly show the interaction of the radio jet with ambient 
gas at systemic velocity. This has not been seen previously 
and indicates that the jet is launched close to the plane of 
the galaxy, pushing the gas encountered on its way. 

• Molecular hydrogen emission arises in extended regions 
along the axis of the galaxy's stellar bar, avoiding the region 
of the bi-cone. The H2 velocities are close to systemic with a 
small rotation component, supporting an origin in the galaxy 
plane. The emitting molecular gas is probably just the "hot 
skin" of a larger reservoir that may have been built up by the 



previously observed HI inflow along the bar of the galaxy, 
and is probably the source of the AGN feeding. 
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